Deterministic Soluble Model of Coarsening 



L. Frachebourg and P. L. Krapivsky 
Center for Polymer Studies and Department of Physics, Boston Umversity, Boston, MA 02215 



> 
O 



o 
o 



We investigate a 3-phase deterministic one-dimensional phase ordering model in which interfaces 
move ballistically and annihilate upon colliding. We determine analytically the autocorrelation 
function A(t). This is done by computing generalized first-passage type probabilities Pn{t) which 
measure the fraction of space crossed by exactly n interfaces during the time interval (0,t), and 
then expressing the autocorrelation function via P„ 's. We further reveal the spatial structure of the 
system by analyzing the domain size distribution. 
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I. INTRODUCTION AND THE MODEL 



We examine phase ordering dynamics in a one- 
dimensional system with three equihbrium states. In our 
model, interfaces between dissimilar domains undergo 
ballistic motion and annihilate upon colliding. The pro- 
cess is thus deterministic although randomness is hidden 
in the initial conditions. Given an appealing simplicity 
of the rules governing the dynamics, it is not surprising 
that this process, and its generalizations, naturally arise 
in different contexts ranging from ballistic annihilation 
to growth processes |6|-||] and dynamics of inter- 
acting populations [|To|-|l2t. Different viewpoints on the 
same model are very useful in that they suggest investi- 
gation of several correlation functions, some of them may 
be clearly interesting from one point of view while could 
hardly be thought from the other point of view. One such 
correlation function, namely the autocorrelation function 
to be determined below, naturally appears in the context 
of population dynamics [p^ ; from other viewpoints, e.g. 
in the original framework of ballistic annihilation ^ it is 
not clear how to define the autocorrelation function. 

We start by describing the two- velocities ballistic anni- 
hilation model and reminding its known basic properties 
The model assumes that interfaces may have two 
different velocities, ±1 without loss of generality, and the 
densities of both populations of interfaces are equal each 
other (otherwise the minority population quickly disap- 
pears). The interfaces are initially randomly distributed 
according to a Poisson distribution. The model exhibits a 
two-length spatial structure, with length scale t[t) ~ \/t 
describing the average distance between neighboring in- 
terfaces moving in the same direction, and the length 
scale C{t) ^ t describing the typical distance between 
neighboring interfaces moving in the opposite directions. 
As we shall see below, however, the growth law for £{t) 
cannot fully characterize the spatial structure - other 
natural measures of the spacing between similar neigh- 
boring interfaces behave differently, e.g.., the rms separa- 



tion £2it) — \J (x^) grows as 

^3/4. We shall argue below 
that all these length scales can be understood as the out- 
come of the competition between the length scale 0(1) 
characterizing initial data and the ballistic length scale 



L{t) - t. 

On the language of phase ordering dynamics, the two- 
velocity ballistic annihilation model may be treated as 
the 3-phase, or 3-color, process with deterministic non- 
conservative dynamics. Indeed, imagine that the one- 
dimensional line is drawn in three colors, say red, green, 
and blue. Let the interface between red and green do- 
mains always moves inside the green one, the interface be- 
tween green and blue domains moves inside the blue one, 
and the interface between blue and red domain moves in- 
side the red one. Then the autocorrelation function A{t) 
is defined as the probability that at a given point and at 
time t the color is identical to the initial color. In the 
dynamics of interacting populations, this model mimics 
a 3-species cyclic food chain [ p^ . 

The rest of this paper is organized as follows. General- 
ized first-passage probabilities are determined in section 
II. Section III contains a calculation of the autocorrela- 
tion function. The domain size distribution is analyzed 
in section IV. The last section V provides a summary and 
an outlook. 



II. GENERALIZED FIRST-PASSAGE 
PROBABILITIES 

Our first goal is to compute P„(i) which measures the 
fraction of space crossed by exactly n interfaces during 
the time interval (0,t). Equivalently, P„(i) is the prob- 
ability that a point has undergone exactly n changes of 
color. Clearly, the color of arbitrary point changes cycli- 
cally with period 3, so the autocorrelation function is 
found from relation 



(1) 



To determine Pn(t), it proves convenient to consider 
an auxiliary one-sided problem with a finite number of 
interfaces on one side of a target point. Namely, imag- 
ine that we have N interfaces to the right of the ori- 
gin (the target point). What is the probability Qn{N) 
that exactly n interfaces will cross the origin? To solve 
for Qn{N), we construct the following discrete random 
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walk: Let Sq = and Si are defined recursively via 
Si = Si-i + Vi, i = 1, . . . , N , where Vi = ±1 is the veloc- 
ity of the i"^ interface. Thus we indeed have a random 
walk {i, Si) starting from the origin, with i being a time- 
like variable and Si a displacement. The crucial point 
is that the number of interfaces which will cross the ori- 
gin is given by the absolute value of the minimum of the 
random walk. Thus we identify Qn{N) with probability 
that an TV-step random walk starting at the origin has a 
minimum at —n. This probability is simply found to be 



Qn{N) ^Qn{N)+Qn+l{N), 



with 



Qn{N) 



m 



2N ^ Af+K ^i ^ ^-" ^1 



(2) 



(3) 



if n and N have the same parity; otherwise, Qn{N) — 0. 

Before returning to the original two-sided problem we 
consider the one-sided problem with infinite number of 
interfaces initially placed to the right of the origin at 
random with density one. During the time interval (0, t) 
interfaces initially located at distances x < t could cross 
the origin. Clearly, the probability Qn{t) that exactly n 
interfaces cross the origin up to time t is 



check of self-consistency we verify that the normalization 
condition, 



(8) 



n=0 



is satisfied. Indeed, Eq. imphes = (X^Q'O^j 

and the latter sum is shown to be equal to one by using 
Eq. d) and identity h{t) ^^Y^^yi I jit) - e* fw). 

Note that P„'s, and especially the first "persistence" 
probability Pq (^), have attracted a considerable recent in- 
terest, see e.g. |15| - |20| . These quantities can be thought 
as first-passage time probabilities in the interacting par- 
ticle systems Given the importance of the first- 
passage type quantities in the classical probability theory 
]l3t one can envision numerous applications of P„ 's in the 
interacting particle systems. However, apart from a few 
findings in the framework of mean-field approach (more 
precisely, for interacting particle systems on a complete 
graph) 1 12 1^ and a limiting analytical solution for the 
ID voter model no exact results are available. The 
model we consider here is an exception in that the com- 
plete analytical solution for P„'s exists, see In 
particular, we have Po{t) — Ql{t) ~ 2(7rt)~^. 

To make results more transparent, it is useful to ex- 
press solutions in the scaling limit 



Qnit) = J2 



m ■ 



Substituting and (||) into (||) yields 



Qn{t) 



[I„{t) + In+l{t)] . 



(4) 



(5) 



where /„ denotes the modified Bessel function of order n. 
If the origin is not crossed by a right moving interface up 
to time t, an interface starting from the origin and mov- 
ing with velocity will survive up to time t/2. Thus 
the surviving probability, S{t), of an interface is given by 



S{t) = Qo{2t) 



-2t 



[Io{2t) + h{2t)]. 



(6) 



First-passage probabilities P„(i) corresponding to the 
two-sided problem are readily expressed via one-sided 
probabilities Qn (t) after realizing that in a configuration 
with n interfaces crossing the origin in the right-sided ver- 
sion, and k interfaces crossing the origin in the left-sided 
version, the total crossing number in the two-sided ver- 
sion is equal to max(fc,rt), as illustrated in Fig. 1. Thus 
we arrive at the relationship 



Pn{t) = 2Q„(t) J2 Qk{t) - 



(7) 



A;=0 



with factor 2 accounting for the fact that smaller number 
k of crossing interfaces can come both from the left and 
right. We have subtracted the last quantity Qn{t)^ which 
has been counted twice in the summation. As a useful 



oo. 



t oo. 



V2t 



= finite. 



(9) 



V nt 

for the one-sided probabilities and 



Making use of the scaling behavior of the modified Bessel 
functions ||||, /„(t) ~ {2'Kt)~'^/^ exp(t - n^/2t), we find 



(10) 



P„(i)=^^^e-^ Erf(^) (11) 
for the two-sided probabilities. 



III. THE AUTOCORRELATION FUNCTION 



The scaling expression of Eq. (|11|) does not allow to 
obtain non trivial long-time behavior of the autocorrela- 
tion function. Indeed, substituting ( [Til ) into (|^) yields 
A{t) ~ 1/3. We should therefore return to exact rela- 
tions (||)-(0). We also extract the trivial A{oo) = 1/3 
factor and consider three autocorrelation functions, 



n=0 



(12) 



describing three possible color outcomes at time t, the 
same (say red) color as initially corresponds to a = 0, 
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Ao{t) = A{t) — 1/3; the "next" blue color corresponds to 
a = 1; and finally the green color corresponds to a = 2. 

All three autocorrelation functions Aa{t) exhibit sim- 
ilar asymptotic behavior; additionally, they are related 
by identity Ao{t) + Ai{t) + ^2(0 = 0. Combining 
and (^) we obtain 



3Ao(i) =3^P3nW-l 

OO 

= [(-^3" - ^3«-l) + (^3" - ^3n+l)] (13) 



n=0 



where F_i = 0. Eq. (Q) allows us to express P„'s via 
QnS. Thus we get 



Psn — Psn-l — Q 



3n "T ^3n-l 



3ra-l 



fc=0 



and 



3n-l 



fc=0 

Substituting Q and (|l|) into yields 

CXD 

3Ao(0 =^(Q3n-Q3n+l)' 

oc 

+ E(Q3„-i+QL-2QL+i) 

OO 3n— 1 

+ 2E(2Q3n - Qsn-l - Qan+l) ^ Qfc- (16) 



ri=0 



fc=0 



In the following calculations we use the exact solution (||), 
the asymptotic relation /„(t) ~ (27ri)~^/^ exp(i — n^/2t), 
and the identity 



2n 

/„_l(i) -/„+l(i) ^—In{t). 



(17) 



The first sum in the right-hand side of Eq. ( |l6| ) behaves 
as 



°c ^-3/2 
/ XQ3n — <33n+l) — „ ^ ■ 



(18) 



n=0 



The second sum in the right-hand side of Eq. (|6|) is 
undergone by treating Q^-iit) — 2Q'^{t) + Qn+iit) as 
the second derivative, d^Q^/dn^ , which is asymptoti- 
cally correct. Using the scaling expression ( p^ for Q„(t), 
this sum is shown to decay as t~'^ in the scaling limit. 
Similarly, the computation of the third line in the right- 
hand side of Eq. (16) is simphfied by the approximation 



Qn-iit) - 2Qn{t) + Qn+i{t) ~ d'^QrJdn^. After some 
algebra, this third term is found to decay as — (2/37r)i~^ 
and thus provides a dominant contribution. The cor- 
responding values for Ai{t) and A2{t) follows from the 
same kind of computation. Thus, wc finally arrive at 
the following asymptotic behavior of the autocorrelation 
bmctions: 

^«W--9^' ^^W-^, Mt)^-^-, (19) 

It is surprising that in the long time limit Aq and A2 
exhibit similar behaviors while the amplitude of the Ai 
has the opposite sign and twice bigger. 

In the general context of coarsening , the autocor- 
relation function is known to decay as >C~^. It has been 
argued that the exponent A satisfies d/2 < \ < d va d 
dimensions Q. Our model implies A{C) ^ Cr^ (A — 1) 
and thus coincides with the upper bound as it happens 
in a few other models, e.^., in the voter model jl9|. Most 
of other studies [|4j also found values of the autocorrela- 
tion exponent satisfying (i/2 < A < d (see, however, Ref. 
p5| reporting the violation of the upper bound for the 
conserved dynamics). 



IV. THE SPATIAL STRUCTURE 

Turn now to the spatial structure formed as the ballis- 
tic annihilation process proceeds. Among several quan- 
tities characterizing the spatial distribution we choose 
the domain size distribution function for which some 
analytical results are already available Let us de- 
note by /i-i (a;,i) the probability density that at time t 

the right nearest neighbor of a -f interface is a — inter- 
face located at distance x apart. Similarly, we introduce 
/i++(x,t) = [1, (a;,i) and [i, \.{x^t\ The Laplace trans- 
form, //(z, €) = Jp°° dx e~^^/x(x, t), of these quantities has 
been computed exactly [0: 



l + J + 2z' 



(z,t) 



S{t)e- 



-2zt 



l + J + 2z' 



fj.+-{z,t) = 



e^^* J^ + 2z(J- 1) 
'Sjrj 1 + J + 2z ' 



(20) 



(21) 



(22) 



where S(t), the probability for the interface to survive 
up to time t, is given by Eq. (^, and 

J=J{z,t) = e-^'*S{t) + 2zf dTe~^''''S{T). (23) 



The solution of Eqs. (|20D-(p2|) has been originally derived 
in an alternative analytical approach to simpler previous 
ones ; this approach of Ref. B has an advantage of 
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being applicable to more difficult ballistic annihilation 
processes like the three velocity ballistic annihilation . 
However, the actual spatial characteristics have not been 
extracted from Eqs. (|ol)-(|2|). 

As a first step, we compute the average length scale 



dx xfi{x, t) 
J^dxfi{x,t) A(0,t) dz 



dfi{z,t) 



(24) 



z=0 



After straightforward calculations we find the average 
size of a domain with boundaries moving in the same 
direction, 



1 + S{t) 



dTS{T)~tS{t) + l 



(25) 



Similarly, we find the average domain size in two other 
situations: 



{x)-+^2t+{x)++, 



(26) 



and 



= 2 



1 - S{t) 
2 

1 + S{t) 



2t- 



Sit) 



drSir) 



dTS{T)-tS{t) + l 



(27) 



Making use of the asymptotic relation S{t) ~ (7r<) -'^Z^, 
we arrive at the following long-time behaviors: 



2t 



2(7r-3)t. (28) 



The latter result is surprising, one might expect that 

(a;)-| grows as (a;)++ while in fact it grows much faster. 

It is instructive to proceed by computing (x")^/" for 
arbitrary positive integer index n. One readily expresses 
via fi{z,t),e.g., (x^) = [/i(0, i)]-i |^^^. 
Any of these quantities can be used to characterize the 
length scale. For domains with dissimilar boundary in- 
terfaces one finds t, implying that 
all these distances are characterized by the single bal- 
listic length scale, C{t) ~ t. In contrast, for simi- 
lar interfaces we get anomalous asymptotic behaviors: 
(x2)i/2 ^ (9^)-i/4i3/4^ generally - t^-^'^-^ 
for integer n. This odd feature indicates the length scale 
characterizing the average separation of the nearest sim- 
ilar moving interfaces, i{t) = ^ is just one 

of the hierarchy of length scales in{t) = (a;")+"- All 
these scales are better thought as effective scales result- 
ing from the competition between the two basic scales 
in the problem, the scale of order one forced by initial 
conditions and the ballistic scale of order t. 

To clarify these results, we compute the inverse Laplace 
transform of Eqs. (|^,|l|,|2|) . We first note that J(z,i) 
can be rewritten as 



/•OC POO 

J{z,t)^2z dTe-^'''S{T)- dTe~'^'^S'{T) 
Jo Jt 

= -z + Vz^ + 2z+ r dre-'^^—^^ 



(29) 



2t 



where we have computed the Laplace transform S{2z) 
and the derivative of S{t). We then expand to find 

fi++{z,t) = a{z) - a{zfh{z,t) + d{zfb{z,tf + ... (30) 
where 



d{z) = z + 1 - ^{z + iy 



1 



and 



6(z,t) 



dre 



(31) 



(32) 



2t 



Performing the inverse Laplace transform of a,{z) and 
6(z, i), we get 

a{x)^^—!^, b{x,t) = ^^^^e{x~2t). (33) 



Combining 

fi++{x,t) = 



OD and (^ 



we finally obtain 



* b 



*2 



»4 



* b 



*3 



(34) 



where f * g — Jq dy f{y)g{x — y) is the convolution of / 
and g, f*^ ^ } * /, f*^ ^f*f* /, etc. 

Noting that a*''{x) = ke~^Ik{x)/x ||2^, the convolu- 
tion in the second term of the right-hand side of Eq. ( |3^ ) 
can be calculated in the long time limit to yield: 



X {l~e-« [/o(0 + 2/i(0+/2(0]} 



(35) 



where = x — 2t. The following terms in Eq. (|3J) give 
corrections for a; > 4i, 6t, . . . 

The only contribution to /x_|_+(a;,i) for a; < 2t is the 
first time- independent term a{x). For large x, it scales 
as a;^'^/^. According to the mapping of Section II, it is 
analogous to the probability that a random walker start- 
ing at the origin first returns back to the origin after 
X steps. A singularity in the second derivative arises 
at a; = 2t, and weaker and weaker singularities appear 
for X being integer multiple of 2t. It should be noted 
that the scale t ^ \/t does not appear in this distribu- 
tion. The only scales appearing are the scale 0(1) char- 
acterizing the time-independent contribution a(a:), and 
the ballistic scale 0{t) characterizing next terms. The 
asymptotic average distance (a;)++ ~ ^At/i: of Eq. ( |25| ) 
is readily obtained by the integration J^^* dx x^++ (x, t) = 

J^'^ dxxa{x), following terms give corrections 0(1). 

Using properties of Laplace transform the dis- 
tributions fj, \-{x,t) and fj, \-{x,t) can be expressed via 

fi++{x,t), 
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/i^ {x,t) = S'(t) /i++(a; — 2<, i) 8(x — 2t) Foundation, and the work of PLK was partially sup- 

n (x + 2t) ported by ARO (grant DAAH04-93-G-0021) and NSF 

^i-+ix,t) = , (36) (grant DMR-9219845). 

thus providing a comprehensive description of the inter- 
faces distribution in this problem. 



V. SUMMARY AND OUTLOOK 

We have shown that the two- velocity ballistic annihila- 
tion process may be thought as the 3-phase deterministic 
model of coarsening. This is one of the simplest models of 
coarsening ever known, and we have derived exact solu- 
tions for the generalized first-passage probabilities P„ (t) , 
and for the autocorrelation function. 

We have revealed a rich spatial structure arising as the 
phase separation process develops. In particular, the mo- 
ments of the domain size distribution, €„(t) = {x")^^, 
exhibit a variety of scales from the time-independent one 
to the scale linearly growing with time: 



.(0 



l-l/2?i 



when n < 1/2, 
when n > 1/2. 



(37) 



We have argued that only the two extreme scales, the 
ballistic one and the scale 0(1) characterizing the initial 
distribution, are important while the others are effective 
in that they arise as the result of competition between 
the extreme scales. The distribution of nearest neigh- 
bors has shown a non-trivial behavior with singularities 
at each x being an integer multiple of 2t. 

Using the mapping on a random walk problem intro- 
duces in Section II, it should be possible to compute 
the two-point equal-time correlation function G{x,t) and 
even the most general two-point correlation function 
C{x,t\0,t') which contains both the equal-time correla- 
tion function G{x,t) = C{x,t\0,t) and the autocorrela- 
tion function A{t) = C(0,t|0, 0). We were able to solve 
for G{x,t) for x > 2t, but the solution is very cumber- 
some so we could not derive clear scaling results. Numeri- 
cal simulations, however, reveal an interesting oscillatory 
behavior of G{x,t). 

Another interesting question concerns the extension of 
the 3-phase deterministic model to higher dimensions. 
It is very simple to define a 3-color cyclic lattice model 
in arbitrary dimension [ p^JT2| . The problem is the sys- 
tem does not exhibit coarsening when d > 2 and instead 
approaches a reactive state with the average number of 
color changes growing linearly with time. However, one 
can hope that a proper higher-dimensional extension still 
exists. 
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